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ABSTRACT 

Context. The detection of GeV photons from gamma- ray bursts (GRBs) has important consequences for the interpretation and mod- 
elling of these most-energetic cosmological explosions. The full exploitation of the high-energy measurements relies, however, on the 
accurate knowledge of the distance to the events. 

Aims. Here we report on the discovery of the afterglow and subsequent redshift determination of GRB 080916C, the first GRB de- 
tected by the Fermi Gamma-Ray Space Telescope with high significance detection of photons at energies >0. 1 GeV. 
Methods. Observations were done with 7-channel imager GROND at the 2.2m MPI/ESO telescope, the SIRIUS instrument at the 
Nagoya-SAAO 1 .4 m telescope in South Africa, and the GMOS instrument at Gemini-S. 

Results. The afterglow photometric redshift of z = 4.35 ± 0.15, based on simultaneous 7-filter observations with the Gamma-Ray 
Optical and Near-infrared Detector (GROND), places GRB 0809 16C among the top 5 % most distant GRBs, and makes it the most 
energetic GRB known to date. The detection of GeV photons from such a distant event is rather surprising. The observed gamma-ray 
variability in the prompt emission together with the redshift suggests a lower limit for the Lorentz factor of the ultra-relativistic ejecta 
of F > 1090. This value rivals any previous measurements of F in GRBs and strengthens the extreme nature of GRB 080916C. 

Key words. Gamma rays: bursts - Techniques: photometric 



1. Introduction 

Long-duration Gamma-Ray Bursts (GRBs) are the high-energy 
signatures of the death of some massive stars and emit the 
bulk of their radiation in the 300-800 keV band. In a small 



number of events, emission up to -lOOMeV has been de- 
tect ed, e.g., with SMM (Harris & Share 1998)| COMPTEL 
(Hoover et al 2005),, EGRET,(Kaneko et al. 2008)| and recently 
with AGILE (Giuliani et al. 2008) These high-energy photons 
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Table 1. Log of the observations. The first ground-based imaging was obtained with the Simultaneous 3-color (JHK) InfraRed 
Imager for Unbiased Su rvey (SIRIUS,|(Na gayama et al. 2003 )[ ) on the Nagoya-SAAO 1.4m telescope (IRSF). GROND, a simulta- 
neous 7-channel imager [(Greiner et al. 2008), mounted at the 2.2m MPI/ESO telescope at La Silla (Chile), started observing about 
30.75 hrs after the GRB. The imaging sequence consisted of a series of sixteen 375 s integrations in the g'r'i'z' channels with gaps 
of about 45 s. In parallel, the JHKs channels were operated with 10 s integrations, separated by 5 s. Late-time imaging was obtained 
with the Gemini-South telescope + Gemini Multi-Object Spectrograph (GMOS-South) on 29 Oct 2008, taking eight 180 s expo- 
sures. Data reduction was done using IRAF routines. Photometric calibration of the GROND g',r', i',z' bands was performed using 
the spectrophotometric standard stars SAlOO-241 and SA97-249, while that of JHKs was done against 2MASS (Tab.|2]i. 



Date/Time 


Telescope/Instrument 


Filter 


Exposure 


Brightness 


(UT in 2008) 






(min) 


(mag)'"* 


Sep 17 02:53-03:43 


IRSF/SIRIUS 


JHKs 


50.0 


21.0±0.5/20.4±0.4/ 20.3+0.5 


Sep 17 07:57-09:39 


MPI/ESO 2.2m/GROND 


g'r'i'z' 


75.0 


>23.6 / 22.81+0.07 / 22.05+0.05 / 21.76+0.05 


Sep 17 07:57-09:39 


MPI/ESO 2.2m/GROND 


JHKs 


60.0 


21.50 + 0.06 / 21.29 + 0.08 / 21.10 + 0.15 


Sep 19 08:04-09:42 


MPI/ESO 2.2m/GROND 


g'r'i'z' 


79.4 


> 23.6 / > 23.8 / 23.47 + 0.13 / > 23.8 


Sep 19 08:04-09:42 


MPI/ESO 2.2m/GROND 


JHKs 


64.0 


> 21.9 /> 21.2 /> 20.5 


Sep 20 08:42-09:42 


MPI/ESO 2.2m/GROND 


g'r'i'z' 


50.0 


> 23.9 / > 24.2 / 23.78 ± 0.16 / > 23.8 


Sep 20 08:42-09:42 


MPI/ESO 2.2m/GROND 


JHKs 


40.0 


> 22.5 /> 21.5 /> 20.6 


Sep 24 07:32-09:31 


MPI/ESO 2.2m/GROND 


g'r'i'z' 


100.1 


> 25.0 / > 24.5 / > 24.3 / > 23.9 


Sep 24 07:32-09:31 


MPI/ESO 2.2m/GROND 


JHKs 


74.2 


> 22.2 /> 21.5/ > 20.8 


Oct 29 07:59-08:31 


Gemini-S/GMOS 


i' 


24.0 


>25.1 



Not corrected for Galactic foreground reddening of E(B-V) = 0.32 mag |(Schlegel et al. 1998)| All magnitudes are given in the AB system. 

Table 2. Local photometric standards within 2' of the GRB. Calibration of the field in JHK was performed using 2MASS stars. 
The magnitudes of the selected 2MASS stars were then transformed into the GROND filter system and finally into AB magnitudes 
using 7(AB) = 7(Vega) + 0.91, //(AB) = //(Vega) + 1.38, K{AB) = /:(Vega) + 1.81 [(Greiner et al. 2008)| Systematic eiTors are 
+0.02 mag for g'r'i'z', and +0.05 mag for JHK^. 



No 


Coordinates (J2000) 


g' 


r' 


i' 




/ 


H 


Ks 


1 


07:59:28.97 -56:38:24.0 


17.59+0.01 


16.82+0.01 


16.50+0.01 


16.27+0.01 


16.24+0.01 


15.98+0.01 


16.36+0.01 


2 


07:59:27.40 -56:40:10.1 


17.18+0.01 


16.18+0.01 


15.79+0.01 


15.52+0.01 


15.37+0.01 


15.00+0.01 


15.38+0.01 


3 


07:59:24.01 -56:37:08.0 


17.11+0.01 


16.31+0.01 


15.98+0.01 


15.72+0.01 


15.67+0.01 


15.40+0.01 


15.77+0.01 


4 


07:59:19.84 -56:39:25.3 


17.90+0.01 


16.92+0.01 


16.54+0.01 


16.28+0.01 


16.15+0.01 


15.84+0.01 


16.20+0.01 


5 


07:59:17.70 -56:37:41.9 


18.11+0.01 


17.09+0.01 


16.71+0.01 


16.45+0.01 


16.36+0.01 


16.02+0.01 


16.44+0.01 



oiTer unique access to the physics of GRBs. Firstly, the shape 
of the spectrum provides direct information about the gamma- 
ray emission mechanism (Pe'er et al. 2007, Giannios 2008). 
Secondly, high-energy photons can place tight constraints on 
the Lorentz factor of the ejecta via the pair-production thresh- 
old. Furthermore, in some cases, the origin of the high- 
energy component differs from that of the low-energy emis- 
sion (e.g., GRB 941017; [(Gonzalez et al. 2003)[ ) or the high- 
energy photons arrive with a significant time delay (e.g., > Ihr 
in GRB 940217; (Hurley et al. 1999) ). The formation of these 
properties is far from understood and can only be addressed with 
an increasing number of bursts with GeV detections. Finally, the 
search for signatures of absorption against the intergalactic UV 
background light using the shape of the high-energy spectrum as 
well as the search for quantum gravity dispersion effects over 
cosmic distances in the hght curve are implications of much 
broader scientific interest [(Abdo et al. 2009)| 

An important prerequisite to the interpretation of the GeV 
component of a burst is the accurate knowledge of the distance. 
Only few of the previously detected GRBs with high-energy 
emission had identified optical afterglows, as the localization ca- 
pabilities of high-energy missions were insufficient to facilitate 
rapid follow-up observations. 

The recently launched Fermi Gamma-Ray Space Telescope 
has the ability to localise high energy events using the 
Large Area Telescope (LAT) and to measure spectra over 



a large energy range in combination with the Gamma-Ray 
Monitor (GBM) (8 keV to 300 GeV). The Swift satellite 
[(Gehrels et al. 20()4)[ can also slew rapidly to LAT locations and 
provide positions with arcsec-accuracy by the detection of the X- 
ray afterglow, facilitating and dramatically enhancing the likeli- 
hood of a distance measurement. 

The bright GRB 0809 16C was det ected by the GBM on 2008 
Sep 16th, at 00:12:45 UT .(Goldstein & van der Horst 2008l[ 
The burst was located in in the field of view of the 
LAT and emission above 100 MeV was quickly localized 
(Tajima et al. 2008l[ Follow-up observations with the Swift X- 
ray telescope (XRT) provided an X-ray afterglow candidate 
|(Kennea 2008), which subsequently led to the discovery of a faint 
optical/ NIR source with G ROND [(Clemens et al. 2008a)[ and 
SIRIUS [(Nagayama 2008)[ Further monitoring in both X-rays 
KStratta et al. 2008)[ and in the optical (Cl emens et al. 2008b)[ es- 
tablished the fading and confirmed the source to be the afterglow 
of GRB 0809 16C. 

GRB 0809 16C w as also detected b y other satelHtes 
in addition to Fermi |(Hurley et al. 2008)| AGILE (MCAL, 
SuperAGILE, and ACS), RHESSI, INTEGRAL (SPI-ACS), 
Konus-Wind, and MESSENGER. The preliminary analysis 
of the GBM integrated spectrum over a duration (T90) of 
66 s results in a best fit Band function ( Band et al. 1993)[ 
with Epeak - 424+24 keV, a low-energy photon index a = - 
0.91+0.02, and a high-energy index /3 = -2.08+0.06, giving 
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a fluence of 1.9x10 erg cm ^ in the 8keV - 30MeV range 
|(van der Horst & Goldstein 2008) The spectral results reported 



by RHESSI and Konus-Wind |(Golentskii et al. 2008)| are in 
broad agreement. 



•3' 



1 



•itRT1.9" 



30" 




Fig. 1. /'-band image of the afterglow of GRB 080916C obtained 
with the 7-channel imager GROND at the 2.2m telescope on 
La Silla / Chile 32hrs after the burst. The circle denotes the 
S'wift/XRT error box. Some of the local standard stars of Tab. 
|2]are labeled. A zoom into the innermost region is shown in the 
bottom-right, with the afterglow (AG) and a galaxy (G) 4" from 
the afterglow labeled. 



The measurements of the high-energy emission from GRB 
0809 16C by the instruments of the Fermi Gamma-Ray Space 
telescope are described in Abdo et al. (2009). Here we report on 
the discovery of the optical/NIR afterglow of GRB 0809 16C, the 
measurement of its redshift, and consequently on the recognition 
of its extreme explosion energy and the large Lorentz factor of 
its relativistic outflow. 



The GRB afterglow 

A comparison of GROND observations from Sep. 17 and 19, 
2008 clearly reveals a fading source inside the S'wift/XRT er- 
ror box (Fig.[B, with coordinates RA (J2000.0) = 07*'59"23^.32, 
Decl. (J2000.0) = -56°38'18'.'0 (a.'5 eiTor). The decay between 
1 .3 to 4 d after the GRB is well described by a single power law 
with ao = 1.40+0.05 (Fig.|2]i, compatible within the errors to the 
X-ray decay slope ax = 1.29+0.09. 

A spectral energy distribution (SED) was constructed us- 
ing the GROND magnitudes from the first night of observations 
(Tab.lTJ. The photometrically calibrated data (Tab.|2]i were cor- 
rected for the foreground galactic reddening of E(B-V) = 0.32 



mag (Schlegel et al. 1998) corresponding to an extinction of Ay 
= 0.98 mag and fit by an intrinsic power law (Fy oc v^^) plus 
three different dust models, as well as without extinction (Tab. 
O. The /' to Ks band data are best fit with a power law slope of 
y6 - 0.38+0.20 and no host-intrinsic extinction. The /' to r' band 
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Fig. 2. X-ray (upper panel) and optical/NIR light curve (lower 
panel) of the GRB 0809 16C afterglow. The solid lines mark the 
best fit power laws to the X-ray and /'-band data (labeled in the 
upper right corner). The power law segments as given by Swift 
[(Stratta et al. 2008)| are shown as dashed lines; they were scaled 
to the /'-band in the lower panel to show that this 3-segment 
power law does not fit the optical data, i.e. the apparent X-ray 
plateau phase is not detected in the optical data. However, the 
X-ray light curve can also be fit with a single power law, with 
only marginally larger as compared to the 3-segment pow- 
erlaw, which results in ax = 1.29+0.09, The optical decay is 
mainly constrained by the three /'-band data points. The JHKs 
observations of SIRIUS and GROND during the first night are 
consistent with this decay. 



measurements deviate significantly and can be best explained 
with a Ly-a break at z = 4.35+0.15 (see Tab. [3] and Fig.[3l and 
the rejected alternative explanations given in the figure caption). 
The redshift values resulting from all the fitted models are com- 
patible, and the redshift error includes already the dependence 
on the error of the photon index as well as foreground-extinction 
correction (Fig.lUi. The g'-band upper limit is consistent with the 
high-redshift result, though it is not deep enough to constrain the 
fit. 

No X-ray measurements of the afterglow are available during 
the time of the initial GROND epoch. However, we can use the 
SIRIUS JHKs -band brightness as well as a back-extrapolation 
of the afterglow decay slope to re-scale the GROND SED to 
the earlier time when XRT measurements are available. The 
Swift/XRT spectrum from 61-102 ks post-burst has been re- 
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Fig. 3. Spectral energy distribution of the afterglow, derived from 
the coaddition of all GROND exposures from the first night (Sep. 
17, 2008). The SED is best fit with a power law of spectral in- 
dex yS = 0.38, no extinction, and Lya absorption at a redshift 
of z = 4.35 + 0.15. Alternatively, the r'-band drop could be at- 
tributed to reddening in the host galaxy, caused either by sub- 
stantial UV absorption or a strong broad absorption feature like 
that at 2175 A |(Kriihler et al. 2008)1 However, the resulting host 
extinction corrected spectral slope of yS S 0.0 would be incom- 
patible with most theoretical models [(Sari et al. 1998)| and with 
the X-ray spectrum. The lack of curvature in the - Ks SED, 
and its extrapolation to the X-ray data, provides additional ar- 
guments against host extinction. Similarly, a spectral break can 
not easily explain the r' - i' color without redshift as the dif- 
ference in power law index would be 2.5, much larger than pre- 
dicted by theor) ](Sari et al. 1998) ' In addition, the steep power 
law would significantly underpredict the observed X-ray fluxes. 
The inset shows that the best-fit GROND power law connects 
without break or offset to the Swift/XRT data, supporting the 
con-ect modelUng of the GROND SED. 

Table 3. Results of the spectral energy distribution fitting with- 
out dust and with dust models from the Milky Way (MW), Large 
Mageflanic Cloud (LMC), Small Mageflanic Cloud (SMC). The 
redshift and (3 errors are at the 2cr confidence level. 



Dust model 


Redshift 




A*"" (mag) 


^red 


none 
MW 


4 as+oi^ 

435+8:1^ 


0.38+j|| 
0.38+ 


0+° 

-QO 


1.04 
1.04 


LMC 
SMC 


4 35+"-" 


34+" " 

"•-'^-0.24 


0+" - 

"•"-0.0 


0.95 
1.04 



ported to be well fit with an absorbed power law spectrum with 
photon index Fx - 2.1^q7 and a column density of A^h - 3.7;^j'j 
X 10^' cm~^ KStratta et al. 2008)| Using the result of no excess 
extinction, we re-fit the X-ray spectrum with the column den- 
sity fixed to the galactic foreground value (A^h = L5 x 10^' 
cm"^), and obtain Fx = L49;^q34, consistent with the slope of 
the GROND SED (note that the spectral index jS is related to the 
photon index F from the X-ray spectral fitting by ySx = Fx - 1). 
The GROND and XRT combined SED is compatible with a sin- 
gle power law over the complete spectral range (see inset of Fig. 

No counterpart or host galaxy was detected seven weeks af- 
ter the burst at the position of the optical/NIR afterglow, with 



S 2.0 




0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



68% 
90% 
99% 




1.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 4.6 4.7 

redshift 



Fig. 4. Contour plot of the foreground Ay (upper panel) and 
spectral index (lower panel) from the spectral energy distribu- 
tion fit against the redshift. This shows that the photo-z deter- 
mination is stable against uncertainties in the correction for the 
galactic foreground extinction, as well as the powerlaw index. 
The 68%, 90% and 99% confidence contours are plotted. The 
shaded area in the upper panel shows the Ay range according to 
|(Schlegel et al. 1998) In the bottom panel the shaded area shows 
the power-law spectral index from the X-ray spectral fits. 



a two-sigma upper limit of /' > 25.1 mag within a 1.1 "-radius 
aperture centered at the afterglow position. This is not surprising 
given the brightness distribution of the known GRB host galaxies 
|(Savaglio et al. 2009)^ the brightness limit places a loose lower 
limit on the redshift of z > 1 . 

We note that the nearest object visible on our images is a 
galaxy at 4" distance to the East. We obtained an optical spec- 
trum of this galaxy beginning at 05:00 UT on 7 November 2008, 
using the Gemini-South telescope + GMOS-S. We obtained two 
spectra of 900 s each with the R400 grating centered at 8000 A, 
a second-order blocking filter in place, and with the slit ori- 
ented to provide simultaneous coverage of the afterglow posi- 
tion. Observations were carried out at relatively high airmass, 
with 1.7" seeing and variable sky background, and consist of 
two spectra of 900 s each. Bias subtraction, flat-fielding, and 
wavelength calibration were performed using the GMOS reduc- 
tion package under the IRAF environment. The wavelength so- 
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Fig. 5. Prompt emission light curve in the 80 ke V-30 MeV energy band as measured with INTEGRAL SPI- ACS at 50 ms resolution. 
The inset shows a zoom of the main peak. Variability on time scale as short as 100 ms is visible - we measure several 6cr flux 
variations relative to the neighbouring data bins on the 100 ms time scale. 



lution was derived using a CuAr lamp spectrum taken immedi- 
ately after the science spectra. Following extraction, the two one- 
dimensional spectra were coadded to increase signal-to-noise, 
achieving a median S/N of 15.8 over the 6000-10,000A wave- 
length range. We clearly detect the continuum of the resolved 
galaxy at wavelengths greater than 6 100 A, placing an upper 
limit on the redshift of this galaxy of z < 4.0, based on the ab- 
sence of Lyman-alpha absorption. This is below the 99% confi- 
dence range for the redshift of GRB 080916C (Fig.|4|, and thus 
this galaxy is not related this burst. 



Discussion and Conclusions 

It is widely believed that long GRBs are produced in the grav- 
itational collapse of a massive star into a neutron star or black 
hole. It has been argued in the past that the observed total energy 
in GRBs requires that the emission is relativistic. GRB 080916C 
is unrivaled by all previous events in this respect: with its ob- 



served fluence in the 8 ke V-30 MeV band (Hurley et al. 2008) 



and redshift of 4.35±0.15, its total isotropic energy release is a 
staggering 6.5x 10^"* erg, corresponding to 4 MqXc^ ! 

The afterglow decay slope is unlikely a post-jet break decay, 
which allows one to assume that the jet break occurs after the 
final XRT observations at 2x10^ s. A break after this time would 
place a lower limit on the jet half-opening angle of 6° 1 +0? 1 (ISM 
surrounding with density of 1 cm"-') or 2°2+0?l (wind medium), 
which in turn implies lower limits for the beaming-corrected en- 



ergy release of (3 .7 ± 0. 1 )x 1 0^^ erg (ISM) or (4.9 + 0. 1 )x 1 0^ ' erg 
(wind). While we cannot distinguish between these two cases, 
we note that the energy for the ISM case is extremely high, ex- 
ceeding the previous record holders GRB 990123 (6 x 10^' erg) 
and 050904 (7 x 10^' erg) by several factors. For the wind 
medium case, a jet break is expected to be very smooth, and 
might have started near the end of our coverage; thus the real 
opening angle may not differ much from our lower limit. On the 



other hand, using our measurements of the decay slope and spec- 
tral shape, the cooling frequency would be still above the X-ray 
band at t ~ 1 d, which is unusual when compared to the majority 
of GRBs: thus the wind interpretation is less likely than the ISM 
interpretation. 

If the emission were non-relativistic, the required photon 
field at the burst location would be optically thick to pair- 
production ("compactness problem"; Ruderman 1975). It has 
been recognized that in addition to the annihilation of pho- 
tons into electron/position pairs, the scattering of photons by 
either the electron or the positron created in the annihilation 
process, contributes to the optical depth of high-energy photons 
[(Lithwick & Sari 200 1)| In fact, this latter limit is in many cases 
more constraining than the pure annihilation limit. 

The most sensitive instrument that detected high-energy 
photons from GRB 0809 16C was the anti-coincidence sys- 
tem (ACS) of d ie spectrometer onboard INTEGRAL (SPI) 
|(Rau et al. 2005)| At its native time resolution SPI-ACS 
recorded at peak more than 1200 counts per 50 ms in the 80keV- 
30 MeV energy range. This allowed the detection of variability 
on time scales as short as 100 ms with high statistical signifi- 
cance (see inset of Fig.|5]). Using eq. 9 of Lithwick & Sari (2001) 
and the photon index of j3 - -2.08 as measured from GBM 
|(van der Horst & Goldstein 2008)| we estimate a lower limit on 
the Lorentz factor of the ejecta of F > 1090. The previously 
highest limit on the Lorentz factor of a GRB with known redshift 
using this method has been F > 410 |(Lithwick & Sari 2001)| for 
GRB 971214 at z = 3.42, for which the additional assumption 
had to be made that the photon spectrum actually extended to 
very high energies. 

Over the last two years, an alternative method to deter- 
mine the initial Lorentz factor has been employed. This in- 
volves using observations of the rising part of optical af- 
terglows to determine when the blast wave has decelerated; 
the corresponding Lorentz factor at the time of the decel- 
eration is expected to be half of the initial Lorentz factor 



6 



J. Greiner et al.: The redshift and afterglow of the extremely energetic gamma-ray burst GRB 080916C 



Fq |(Sari & Piran 1999)| T his method provided Tq x 400 for TajimaH., Bregeon J., Chiang J., et al. 2008, GCN #8246 

GRB 060418 and 060607 |(Molinari et al. 2007)1 To = 160 for vanderHorst A., Goldstein A., 2008,gcn #8278 
GRB 070802 KKruhler etal. 2008)| Tp ^ 120 for GRB 080129 
[(Greiner et al. 2009)| and Tq = 200 for GRB 071031 (Ki-uhleret 
al. 2009). Yet another method is based on the evolution of the 
thermal emission component in the prompt emission of GRBs 
|(Pe'er et al. 2007) } and also provides similarly low values of F. It 
is interesting to note that our lower limit on F for GRB 0809 16C 
is substantially larger than values determined by other meth- 
ods. Whether or not this is related to the GeV emission in 
GRB 0809 16C remains to be seen. 

The high signal-to-noise ratio of the SPl-ACS data is also 
ideal for estimating the variability of the light curve, a quan- 
tity that has been shown to correlate with the isotropic equiva- 
lent peak luminosity, Liso.peak- Following the method described 
in Li & Paczyiiski (2006) and using a smoothing time scale 
of 13.7 s, we derived a variability index of V = -2.26 and 
a resulting Liso.peak = (1-23 + 0.32) x lO^^erg s"' (80keV- 
30MeV). Using the observed 256 ms peak flux from Konus 
[(Golentskii et al. 2008)| we derive Liso.peak - 2x 10^^ erg s"' 
(20keV-10MeV). 

For the future, the synergy of the detection of GRBs with 
GeV emission coupled with the ability to localise and determine 
redshifts for these events will be extremely interesting as both F- 
determination methods can be applied, providing a consistency 
check of our picture of the GRB and afterglow emission process. 
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